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Abstract 

In this work we revisit the concept of chemical potential u in both classical and quantum gases 
from a perspective of Equilibrium Statistical Mechanics (ESM). Two new results regarding the 
equation of state fi = u(n,T), where n is the particle density and T the absolute temperature, 
are given for the classical interacting gas and for the weakly-interacting quantum Bose gas. In 
order to make this review self-contained and adequate for a general reader we provide all the basic 
elements in an advanced-undergraduate or graduate statistical mechanics course required to follow 
all the calculations. We start by presenting a calculation of fJ,(n, T) for the classical ideal gas in 
the canonical ensemble. After this, we consider the interactions between particles and compute the 
effects of them on /J,(n, T) for the van der Waals gas. For quantum gases we present an alternative 
approach to calculate the Bose-Einstein (BE) and Fermi-Dirac (FD) statistics. We show that 
this scheme can be straightforwardly generalized to determine what we have called Intermediate 
Quantum Statistics (IQS) which deal with ideal quantum systems where a single-particle energy 
can be occupied by at most j particles with ^ j ^ N with N the total number of particles. In the 
final part we address general considerations that underlie the theory of weakly interacting quantum 
gases. In the case of the weakly interacting Bose gas, we focus our attention to the equation of 
state u = fj,(n,T) in the Hartree-Fock mean- field approximation (HF) and the implications of 
such results in the elucidation of the order of the phase transitions involved in the BEC phase for 
non- ideal Bose gases. 
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I. INTRODUCTION 



Chemical potential has proven to be a subtle concept in thermodynamics and statistical 
mechanics since its appearance in the classical works of J.W. Gibbs 1 . Unlike thermodynamic 
concepts such as temperature T, internal energy E or even entropy S, chemical potential 
/j, has acquired, justified or not, a reputation as a concept not easy to grasp even for the 
experienced physicist. Gibbs introduced chemical potential within the context of an extense 
and detailed exposition on the foundations of what is now called statistical mechanics. In 
his exposition he considers how to construct an ensemble of systems which can exchange 
particles with the surroundings. In such description, /i appears as a constant required 
to provide a necessary closure to the corresponding set of equations 1 . A fundamental 
connection with thermodynamics is thus achieved by observing that the until-then unknown 
constant fi is indeed related, through first derivatives, to standard thermodynamic functions 
like the Helmholtz free energy F = E — TS or the Gibbs thermodynamic potential 
G = F + pV. In fact, fi appeared as a conjugate thermodynamic variable to the number 
iV of particles in the same sense as pressure p is a conjugate variable to volume V. The 
procedure outlined above and described with detailed elegance by J.W. Gibbs defines the 
essence of the chemical potential in statistical mechanics and thermodynamics. 

The link provided by Gibbs to define chemical potential in terms of thermodynamic 
variables is certainly a master piece, however, a direct physical interpretation might still 
be elusive. Consider for example, two of the most used definitions of \i in equilibrium 
thermodynamics 2 , 3 , i.e., 



where V is the system volume. As can be readily seen, the first definition in terms of F's 
derivative implies that we can obtain [i as a measure of the change of F with respect to the 
number iV of particles at constant volume and temperature. It is straightforward to imagine 
a closed box of volume V where we can add or subtract particles and observe changes in 
the free energy of the system. However, depicting such situation when both volume V and 
temperature T are kept fixed may require a higher degree of physical intuition recalling 
that any particle added to the system will provide some additional amount of energy either 
in the form of potential or kinetic energy. Let us consider the second term in Eq (1). This 
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thermodynamic definition suggests that fi can be measured as the change of internal energy 
E with respect to the number iV of particles but this time keeping constant entropy S and 
volume V. What exactly must be understood by adding particles keeping entropy constant? 
Recall that each particle added to the system brings an increase in the number of configura- 
tions available to the overall system and therefore an increase of entropy would be our first 
intuitive expectation. In equilibrium statistical mechanics (ESM) the procedure to present 
fi is very much based on the approach suggested by Gibbs in his classic works. The central 
idea behind ESM is a many-variable minimization process in order to obtain a distribution 
function {n q } corresponding to an extremal, minimum or maximum, of the thermody- 
namic variables F or S respectively. In this context fi appears as a mathematical auxiliary 
quantity identified with a Lagrange multiplier that minimizes/maximizes a physical quantity. 

In this work we present a discussion on the concept of chemical potential from a 
perspective of ESM and how it emerges from physical considerations in both classical and 
quantum gases. Our main focus is to present to undergraduate and graduate students a 
self-contained review on the basic elements that give rise to an understanding of fi. In 
order to achieve this goal we shall proceed as follows. Section II carries out a detailed 
calculation of /i in the case of an ideal classical gas. Using a method proposed by van 
Kampen 4 we include the effect of interactions and calculate /i for the van der Waals gas. 
Section III deals with ideal quantum gases. We introduce calculations by giving a general 
discussion on the temperature scale at which quantum effects are expected to contribute 
significantly. A general description, based on a simple, but novel method to compute the 
average number of particles (n^) that occupy the single-particle energy level for boson 
and fermions is introduced. We also provide a formal calculation for (nk)j when the energy 
levels can be occupied at most by j particles, where ^ j ^ N. We call the resulting 
statistics the Intermediate Quantum Statistics (IQS) of order j, which generalizes the BE 
and FD statistics which are obtained for j — > oo and j = 1, respectively. Finally in Section 
IV we go a step further to consider the behavior of /i as a function of particle density n 
and temperature T for the weakly interacting quantum Bose and Fermi gas. The former 
system has been under intense research lately since it is the standard theoretical model to 
describe Bose-Einstein Condensation (BEC) in ultracold alkali atoms 5 , and as we outlined 
here, the knowledge of fi = fi(n,T) is of fundamental importance since it contains valuable 



3 



information on the nature of the phase transition involved. 



Our intention is not to provide an exhaustive treatment of the chemical potential in 
ESM, instead, our contribution intends to integrate previous well known results within a 
physically intuitive framework, and at the same time to provide some new results that 
might be interesting to the reader that complement and enhance a broader view of the 
subject. We kindly invite to the interested reader to study several excellent textbooks 6-9 
and reviews 10 , 11 , 12-14 that have been written on the subject in the recent past. 



II. CHEMICAL POTENTIAL I: THE CLASSICAL GAS 

In order to discuss chemical potential for the ideal classical gas we shall address some basic 
considerations and implications that \i must satisfies according to general thermodynamic 
principles. Although we present them at this point, its validity goes beyond the classical 
ideal gas. Fundamental postulate in equilibrium thermodynamics 2 assures that for a given 
system there is a function called the entropy S defined only for equilibrium states which 
depends on volume V, internal energy E and number of particles N, i.e, S = S(E,V, N). 
Thus, an infinitesimal change dS between two equilibrium states can be written as 

dS= m dE+ m *, + m m. (2 > 

\dEj VN \dVj EN \dNj VE 

Using the first law of thermodynamics dE = TdS — pdV we can relate the partial derivatives 
that appear in Eq (2) with standard thermodynamic variables temperature T and pressure 
p. A simple inspection points out that 

(-) =~ 

(4) 

Such identification suggests that we must add to the First Law a suitable thermodynamic 
variable that will play the role of a conjugate variable to the number of particles N and that 
will allow the connection to (dS/dN) VE , just like T is conjugated to the entropy S and p 
is to the volume V. Thus, if we allow the exchange of particles, we can write the First Law 
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as dE = Tds — pdV + fidN and hence 

Eq. (5) provides additional information on the nature of chemical potential complementing 
Eqs. (1). This tells us that fi is a negative quantity if entropy increases with the number 
of particles by keeping energy E and volume V constant. Though it is intuitive that S 
increases as N increases, it is not the case under the restrictions of E and V constant. On 
the other hand, Eq. (5) also admits the possibility that fi > 0, however, as we show below 
for the ideal Fermi gas and the weakly interacting Bose gas in sections III- IV, respectively, 
this is true only as a result of quantum effects. 



A. The classical ideal gas 

To determine \x as a function of (E, V, N) we shall make use of the fundamental equation 
S = S(E,V,N) and Eq. (5). ESM ensures that the macroscopic variable entropy S is 
related to a microscopic quantity Q(E, V, N) which represents the number of microstates 
available to the system consistent with the macroscopic restrictions of constant E, V and 
N. Such connection is given by S = ks In Q(E, V, N) where ks is the Bolztmann's constant. 
With this considerations, Q is given by 

Q(E,V,N) = -L-J...Js{E - F)d 3 ri d 3 p?d 3 r 3 ivrf 3 P3iv (6) 

where H = Yli=iPl/^ m * s the Hamiltonian for a system of N free particles and 1/N\ corre- 
sponds to the Gibb's correction factor. Since H is r-independent and spherically symmetric 
respect to momentum coordinates Pi we can write Eq. (6) as 

v n 2tt 3N / 2 r°° 

w) - m^TWrn x /„ dpp S(E - p2/2m) (7) 

where the change of variable P = Y^=i Pi ^ as been made and the hyper- volume element in 
3./V dimensions with coordinates pi given by dli = 2n 3N / 2 P 3N ~ 1 dP, with P defined as before, 
has been used. Then, 

1 yN E13JV/2-1 
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In the limit where N 3> 1, Eq. (8) can be written as 

Thus, entropy 5* can be readily calculated. After using Stirling's approximation one obtains 

mE 



S A, ; .Y>^ + ^ln 



3nh 2 N 



Substitution of (10) into (5) leads to the well known result for the chemical potential for the 
ideal classical gas 

V ( mksT^ 3 



= -k B T In ^ v 27r/ _ 2 
where the relation = fiVfeeT has been used. 



(11) 



Eq. (11) has an interesting interpretation in terms of the average distance between 
particles I = (V/N) 1 ^ 3 and the thermal-wavelength At = hj 1 \f2nmkBT (see section III for 
a larger discussion). As it can be seen from Eq. (11), fi = — fc^Tln [/ 3 /A^] from which 
a physical interpretation can be easily harnessed. The sign of the chemical potential is 
determined then by the ratio //At- In the high-temperature limit, when the quantum effects 
are small and the wave nature of particles is negligible in comparison to I, i.e, At *C I, \i 
is negative and the system can be regarded as formed of idealized punctual particles that 
can be distinguished, in principle, one from each other. This picture corresponds to the 
ideal classical gas. This interpretation opens up the possibility that in the quantum regime, 
At ~ I, li could acquire positive values. 

We can gain additional information if we consider the discrete version of Eq (5), namely 

(M = -T(AS) Ey (12) 

where 

Notice that Eq(13) gives the sign of the chemical potential when one particle is exactly 
added to the system keeping E and V constant. Substitution of Eq (8) into Eq (13) yields 

tt(E,V,N + l) V (§iV-l)!/ m \V 2 



Sl(E,V,N) (JV + 1)(|JV + I) 



6 



For iV > 1, 



Q(E, V,N + 1) ^2V 
Q(E, V, N) ~ 3iV 




mE 



) 



3/2 



37ch 2 N 



(15) 



where Stirling's approximation has been used. In this representation, // goes essentially 
as the logarithm of the ratio between the energy per particle E/N and the energy e = 
ft, 2 /2m(V/N) 2 / 3 of a quantum particle confined in a box of volume V. The condition E/N ^> 
e guarantees the classical character of the system assigning a negative value to the chemical 
potential. 

B. The effects of interactions 

It has been shown in previous section that chemical potential for the ideal classical gas is a 
negative quantity for the whole temperature region where quantum effects can be neglected. 
In order to enhance our intuition on the nature of chemical potential we shall address the 
calculation of \x in the case of a classical gas with pairwise interactions between particles. 
For a system of N particles, the total partition function can be written as 6-8 



is known as the configurational integral. In Eq. (17) Vi j = u(|r$ — r 3 -|) is the interaction 
energy between the i-th and j-th particles and (3 = (fc^T) -1 as usual. 

A simple method to evaluate Qn has been given by van Kampen in Ref. 4 . In such work, it 
is suggested that the average of e~P Vl > 2 e~P Vl < 3 ■ ■ .e~P Vlf - 1 > N over all possible configurations of 
particle's positions can be identified exactly as the ratio Qn/V N ■ Then, the configurational 
partition function Q N can be expressed as 4 




(16) 



where 




(17) 




(18) 



where the coefficients are given by 




(19) 
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and the sum is taken over all irreducible terms that involve A;-particle position coordinates 
(see Appendix for more details). The total partition function Zn is then given by 



V N (mk B T 



3N/2 



exp 



k=l 



Bi 



k + 1 



(20) 



and from this, the Helmholtz free energy F by 

V fmk B T\ V2 

N 1 27T^ 2 ) 



F = -« B T>n^=) j - Nk B T 
The chemical potential fi can be obtained readily as 

\i> Hideal 



k=l 



Bi 



V k+l 



■ (21) 



k=l 



Bh. 



(22) 



Eq. (22) gives \x for the classical interacting gas as a series of powers in the particle density 
(N/V) k . The physical implications are clear, interactions shift the value of the chemical 
potential from the ideal case. If there is no interactions at all, then B/, = for all k and 
A* = ^ideai- In spite of the generality of expression (22), in practice, calculation of B k for 
k > 2 is rather cumbersome. However, for enough dilute systems, i.e, N/V 1, we may 
consider only the first term of Eq. (22) as a valid approximation. Thus, at first order in N/V 
we have /x = fiideai — kBT{N/V)B\, where B\ depends on the specific interatomic potential 
between particles. In order to obtain quantitative results about the effects of interactions 
on the chemical potential we consider the van der Waals g specific example. 

Let us consider, for simplicity, the commonly-used pairwise interaction potential 

for r < d 
-vo(d/r) 6 for r > d, 

that approximates the semi-empirical Lennard- Jones potential v(r) = t> [(d/r) 12 — 2(d/r) 6 ] , 
vq is the minimum interaction energy between a pair of particles and d their separation at 
which such energy takes place. For this interaction model, B\ can be evaluated exactly 
as follows. By taking advantage of the spherical symmetry of the problem we can write 



v(r) 



oo 



(23) 



Bi = lie 



3v{t) _ 



l)dr = 4tt/ ( 



00 ^2 



re 



3v(r) 



l)dr. Then, by splitting the last integral into 



one integral from to d plus a second one from d to oo and using the fact that v(r) — > oo 
for < r < d, we get 



Bt = 4tt 



,Pv {d/rf 



d 3 

\)r 2 dr — 



(24) 



8 



The integral in Eq. (24) can be evaluated directly by using the Taylor series of the expo- 
nential function. After integrating term by term we have 



4 , 
B l = -ixd 3 



(fob)" 1 
^ (2n - l)n! 

n=l v ' 



(25) 



It is possible to go a step further in order to write Eq. (25) in terms of elementary functions, 
certainly, the infinite sum can be written as Y^mLi x ^ V $ n ~ = 1 ~ + (7rx) 1 / 2 Erfi(a; 1 / 2 ), 
where Erfi(^) = —i Erf(iz) denotes the imaginary error function. A simple expression for 
the correction of \iideai due to interactions, defined as A/i = jjl — fiideah can De obtained 
for temperatures such that k B T ^> vq, since only the first term in the series expansion in 
expression (25) is needed, with this approximations and recalling that I = (V/N) 1 / 3 we have 



^ aW rlij l 1 -^j >0 (26) 

in agreement with Monte Carlo calculations obtained previously by other authors 15 . For even 
higher temperatures, vo/ksT w 0, and then it is the hardcore repulsion of the inter-particle 
interaction (23) what governs the dynamics of the gas. In this limit the system corresponds 
to a hard-sphere gas thus giving A/i = kBT^n(d/l) 315 . For temperatures smaller than vo/ks, 
Afi becomes negative (see Fig. 1), but this should not be considered correct since at such 
temperatures we are out of the classical regime and quantum corrections must be taken into 
account. In terms of the parameters a and b of the standard van der Waals equation of state 

N 2 \ 

P + ^ a ) (V-Nb) = Nk B T, (27) 

the chemical potential for the van der Waals gas can be written as 

N 

V = l^ideai - 2— (a - k B Tb) (28) 

where a = v b = Vo^nd 3 . Table I presents some standard values for the a and b values 
for different gases 17 . The interested reader may find useful to see how the calculation just 
presented works, by using other interaction potentials Vij between particles. 

A relation of A/i with the work W(r) required to bring an additional particle to the 
system from infinity to position r, has been shown by Widom 16 as 

exp (-Afi/k B T) = (exp(-W(r)/k B T)} } (29) 
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TABLE I: Values of the van der Waals parameters a and b for some substances are given. With 
these values the ratio vq/UbTr is computed, where Tr denotes the room temperature. 



v(r) 







d 



Substance 


a 


b 


Vn/kn Tr 


Helium 


0.0346 


0.0238 


0.0603 


Neon 


0.208 


0.0167 




Hydrogen 


0.2452 


0.0265 


0.384 


Oxygen 


1.382 


0.0319 


1.796 


Water 


5.537 


0.0305 


7.527 



Hard sphere limit 
''(high temperature regime) 




(low temperature regime)\ 



v / k B T 



FIG. 1: Left panel shows the model potential given by expression (23) as an approximation to 
the more realistic Lennard-Jones potential. In the right panel we present the change in chemical 
potential (26) as function of the ratio of the energy that characterizes the interacting potential vq 
to the thermal energy. 

where (..) denotes the canonical-ensemble average. On the other hand, it seems intuitive to 
expect ^(r) to be larger for a gas with repulsive interactions than for the ideal gas, thus, 
by using Widom's equivalence Eq (29) we may conclude that repulsive interactions yields 
A/i > 0. 
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III. CHEMICAL POTENTIAL II: QUANTUM IDEAL GASES 



At low temperatures gases depart from their classical behavior since quantum properties 
of matter emerge. One of the main properties washed away in the quantum regime is distin- 
guishability. In the classical picture, we can in principle label and tag any of the particles, 
but no longer in the quantum regime 18 . This property has profound consequences in the 
number of different microstates available to the system. In general, classical systems will 
have more microstates since permutations among particles result in different configurations 
due to distinguishability. Quantum systems on the other hand display a smaller number of 
different configurations. In addition to indistinguishability, quantum gases exhibit another 
remarkable property. L. de Broglie suggested that any material particle with mass m and 
velocity v should have a corresponding wavelength A given by 

h 2-nfi , , 

A = - = , 30 

p mv 

where h is the Planck's constant and p the momentum of the particle. Given the fact that a 
particle with kinetic energy mv 2 /2 has an associated temperature T, it is possible to write 
down an expression for a thermal de Broglie wavelength Ay as 

Xt = fAttt - (31) 

Eq. (31) establishes indeed a criterion that determines whether the nature of a system of 
particles can be considered as classical or quantum. Basically, the wavelength At serves as 
a length scale over which quantum effects appear. For high temperatures Xt — > and then 
the particles can be visualized as classical point-like particles with a definite momentum and 
position. However, as temperature is lowered, Xt starts to increase is a smooth way. There 
exist then a characteristic temperature T*, such that the wavelength of particles is of the 
same order of magnitude as the average distance I between any two particles (see Fig. 2), 
i.e., 

I ~ A*. (32) 

At this temperature T*, the system enters into the so called, degeneracy regime. In such 
conditions the wave-like properties of matter drive the phenomenology of the system. Eq. 
(32) is much more than a qualitative description, assembled together with Eqs. (30) and 
(31), provide the correct order of magnitude for the critical temperature of condensa- 
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tion T c in ultracold alkali gases used in current experiments of Bose-Einstein Condensation 18 . 

Indistinguishability of particles in the quantum regime requires the A-particle wave- 
function ^(ri, . . . , r/v) of the system satisfies certain symmetry properties. These symmetry 
requirements for the wave function of the A-particle system implies the existence of two 
fundamental classes of quantum systems 19 . A system for which the total wave function is 
symmetric with respect to the exchange on the positions of any two particles, i.e. 



and other system where the wave function is anti-symmetric with respect to this action, i.e, 



The first case corresponds to a system formed by particles called bosons while the second to a 
system formed by fermions. In addition, the expression (34) serve as the basis for the Pauli's 
exclusion principle: no two identical fermions can ocupy one and the same quantum state. 
Both systems, Bose and Fermi gas, exhibit completely different macroscopic properties as 
we shall show below. 

The same symmetry considerations on the wave function has also a direct consequence 
on the spin of the particles involved 19 . It can be shown that for a quantum system with 
a symmetric wave function \1/ (bosons), in the sense described above, the particle's spin s 
can only have integer values, i.e., s = 0, 1,2, . . .. For system with an antisymmetric wave 
function (fermions), particles can only have a spin with positive semi-integer values, that is, 
s = 1/2,3/2,.... Such difference in spin values shall manifest in larger differences in their 
macroscopic dynamics. 

A. The Bose-Einstein Distribution 

In this Section we shall review the thermodynamic consequences of this symmetry 
condition for the chemical potential fi in an ideal Bose gas. If we denote by to the 
number of bosons that populate a particular energy level k, then the symmetry principle 
presented above implies njt = 0, 1, 2, . . . oo. A direct implication of this fact is that at zero 
temperature, a system of Bose particles will have a macroscopic occupation of the lowest 




(33) 




(34) 
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FIG. 2: In the high temperature limit a), quantum statistical correlations measured by the thermal 
wavelength A are much smaller than the average separation between atoms I. b) As temperature 
is lowered, quantum effects start becoming important when relation (32) holds. At even lower 
temperatures c) indistinguishability is dominant. 

single-particle energy level. This was first recognized and discussed by Bose and Einstein 
in 1925 and gave rise to a large interest that culminated in the experimental realization of 
the first Bose-Einstein Condensate in 1995 by Ketterle et al 20 and Weiman and Cornell 21 
with atoms of 26 Na and 87 Rb, respectively. The phase known as Bose-Einstein Condensate 
(BEC) corresponds to the state where the number of particles in the lowest-energy level no 
is of the order of the total number N of particles. 

We start our discussion by noticing that two distinct situations must be considered: 
a case in which the number of bosons A" is conserved at all temperatures and the case 
in which is not. A well established principle in physics tells us that the number of the 
fermions (baryons and leptons) involved in any physical process must be conserved 22 . Thus 
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it is expected that composite atomic bosons, such as the trapped atoms used in the BEC 
experiments also satisfies this principle. There are, however, many physical situations where 
the number of bosons in the system is intrinsically not conserved. Collective phenomena 
that emerge from interacting ordinary matter give rise to bosonic pseudo-particles ("quasi- 
particles" or simply "excitations"), that are created from the system's ground state by 
simply raising the system's temperature for instance. At equilibrium, quasi-particles are 
created and annihilated with a very short life-time due to thermal fluctuations. This is the 
main reason why they are not conserved. Interestingly, it is possible that systems composed 
of quasi-particles can exhibit Bose-Einstein Condensation. Indeed, experimental evidence 
of BEC of spin-excitations (magnon gas) 25 and of particle- hole coupled to photons in an 
optical cavity (exciton-polaritons) 24 ' 26-28 has been reported. The possibility of BEC in this 
systems relies on the fact that the quasi-particle life-time is much larger than the relaxation- 
time, thus the system can be described by statistical mechanics where particle number is 
approximately conserved. 

A important system composed by bosonic particles that are not conserved is electromag- 
netic radiation in a cavity at thermodynamic equilibrium (black body radiation). From a 
quantum perspective, electromagnetic radiation can be considered as an ideal gas composed 
by photons that obey the rules of symmetry of the wave function As in the case of 
quasi-particles, the number N of photons is not conserved since they are continuously being 
absorbed and emitted by the cavity's walls. There is, however, a very recent work under 
progress where an experimental setup has been achieved in order to produce an electromag- 
netic radiation system with conserving number of particles 30 . In the next paragraphs we 
shall address the behavior of the chemical potential for a system of bosons where the number 
of particles is not conserved, as it occurs in the case of electromagnetic radiation confined 
in a cavity. In order to explore the capabilities and strengths of the ESM approach we shall 
partially follow the procedure suggested by Reif 7 . We will find this very instructive when 
we extend this method to deal with quantum system with fractional statistics. 

1. The photon gas 

Let us consider a quantum ideal system composed of N bosons that can be distributed 
along a set of {e^} energy levels. Due to the fact that we are dealing with bosons, each 
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energy level e k might be populated with n k particles with n k = 0,1,2,.. .00. The average 
number of particles (rife) that occupies the single-particle energy level e k as a function of 
temperature is given by the following general prescription 



(n k ) = = , (35) 

where two types of sums are defined at different levels. For a particular macrostate M, there 
is a distribution of n\ particles in the energy level E\, ri2 particles in the energy level e 2 and 
so on. The total energy E of this particular macrostate is given by 

E = Y t n s e j (36) 
3 

where j = 1, 2, . . . , 00 if there is any possibility that the system may be excited to higher 
energy levels as desired. In practice the sum defined in Eq (36) ends at some finite point 
beyond of which there is no possibility that any particle can occupy higher energy levels. 
However this sum is defined only for a particular macrostate. To take into account all 
different possible macrostates M consistent with the same total energy E one must to define 
a second sum on a higher layer. The outer sum Yl{n k } takes into account this fact and it 
must be performed over all possible macrostates M availabe to the system. Thus in spite 
its apparent simplicity Eq. (35) is in fact a sum over the possible distributions of particles 
in the all possible energy levels. In order words, Eq (35) comprises a sum over the different 
sets {ni,n 2 ,n 3 , . . .} compatible with the restriction that for any of these distribution sets 
it must occur that the total energy E of the system must be given by Eq (36). 

The key element to perform the sum defined (35) is to define the partition function Z of 

the system as 

Z{V,T) = Y,e-^ n ^ (37) 

{n k } 

and then split the sum in Z into the state k and the remaining ones. This can be written as 

Z(V T) = e -/3(niei+n 2 e2+— +n fe e fe + ) 

{n k } 

— y ^ e ~Pn k e k g — 0(71161+71262+— ) ^gg^ 

n k n q ^n k 
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where the last sum explicitly excludes the term n^k which has been separated and brought 
up to the front of the sum. Then the quantity (nk) can be written as 

E„, p -Pn k ek V p -/3(n 1 e 1 +n 2 e 2 +...) 

\n k ) - -/3n k e k y e -j9(mei+n2ea+...) ^ ' 

£-/n k Z-^n q ^n k 

where the numerator has only one sum that depends on n k . Up to this point these results 
are completely general. They do not depend on the particular values that the numbers rij 
can assume. 



For bosons the numbers rij can assume any of the values 0, 1, 2, ... N since there is no 
restriction on the energy level occupancy. In the case of electromagnetic radiation and other 
quasi-particle systems there is no constraint on the total number of particles N and thus the 
sums over the states n q ^ n& in Eq (39) are identical and can be canceled out. This yields 
to a simpler expression for the average occupancy, i.e, 

<»*> = x»' ' • (40) 



*—>n k 



In order to carry out this calculation let us define the quantity z(V,T) as 

z(V,T) = J2 e ^ nk€k ( 41 ) 

n k 

It is straightforward to see that in terms of z(V,T) the average number (rife) of occupancy 
can be written as 

Since the numbers can adopt any possible value between zero and infinity, the sum 
expressed in z(V,T) is indeed a geometric series which can be readily calculated as 

oo 1 

,(V,T)=^e-^ = r -^, (43) 

n k =0 

and thus 

= i^TI- ^ 
represents the average number of particles occupying the energy level e k in the case of 
non conserving total number iV of particles. We can apply this result to the photon 
gas. In such case the energy of a photon with wave vector k in a given polarization 

1(3 



is determined by = hu^ = hc\k\. Substitution of this in in Eq. (44) enables us 
to recover the well known Planck distribution which gives the average distribution of 
photons in the uy_ mode. It is of paramount relevance to realize that in this case /i does 
not appear explicitly during the calculation. We might say that chemical potential in 
this case is zero. However, a stronger assertion can be made. Within canonical ensem- 
bles in ESM there is no need of chemical potential if the number of particles is not conserved. 

2. The ideal Bose gas 

Let us now implement the above procedure for a system of bosons where the total number 
N particles is conserved. In this category falls a vast set of system made of actual massive 
boson particles. The only requisite is that the total number of particles N must be kept 
fixed. In order to calculate the average number (rife) in this case we must return to general 
expression for (n^) given in Eq (39). For the case where the number of particles is allowed 
to fluctuate we noticed that the sums ,„ e -^(niei+n 2 e 2 +...) j n ^ e numerator and denom- 
inator were indeed the same and thus they can be canceled out. For systems where there is 
a restriction on the value of N this is not longer the case. To understand this we just must 
realize that the restriction 



constrains the sums given in Eq (39), to be performed over the remaining particles once 
the energy level k has been occupied. This is, if from N particles one is occupying the 
energy level k, then the remaining N — 1 particles have to be distributed necessarily over 
the energy levels q with q ^ k. This simple and powerful idea is the core of the calculation 
presented by Reif 7 . We shall use and extended this idea to calculate averaged occupancies 
(rifc) beyond the BE and FD statistics. 

For a system composed of bosons there is no restriction on the total number of particles 
that can occupy a single-particle energy level e&. Consider for instance, the energy level e q 
being occupied by one particle. Then, Eq. (45) implies than the remaining N — 1 particles 
must distribute themselves into the energy levels k with k ^ q. In terms of Eq. (39) this 




(45) 



k 
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also implies that the sum 

e - / 8(mei+n a e2+...) 

must be carried out not over the all the possible values of n 1; n 2 , etc but only over those 
that satisfy the fact that the total amount of particles available is now N — 1 for this case 
and not N as it was at the beginning. Let us denote this new sum as 

Z'(N-l)= e-P {niei+n2e2+ -\ (47) 

where the prime / denotes the fact that the sum must be performed over all energy levels 
different form q. The quantity iV — 1 in parentheses indicates too that this sum is carried 
out over N — 1 particles. 

With these elements we are in position to calculate (rik) for a system composed of N 
bosons with iV fixed. For each value adopted by nk, the sum in Eq(47) has to be carried 
out over the remaining particles. This yields to 

N 

(n k ) = ^ ■ (48) 

Y,e-^Z\N-j) 

3=0 

The evaluation of Eq (48) requires to compute Z'(j) from j = 1 to N, which makes the 
calculation rather cumbersome, see Ref. 31 and 7 for details. We shall detour this difficulty by 
taking an alternative approach to that followed by Reif. Let us Z S (N) from the numerator 
and denominator of Eq (48) to obtain 

N 



je 



Z'(N-j) 



Z'(N) 

(nk) = ^ (49) 

Z'(N-j) 



°~ ' ; ' ; 

3=0 



6 Z'(N) 



Note that the ratio Z'(N — j)/Z'(N) can be written as the product of the ratios of partition 
functions that differ only in one particle, i.e., 

Z'(N-j) Z'(N-l) Z'(N-2) Z'(N-j) 

Z'(N) Z'(N) Z'(N-l) Z'(N-j + l)' { ' 
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We use now that the finite change of the Helmholtz free energy AF, when exactly just one 
particle is added to an iV-particle system corresponding to a chemical potential 

^ = k ° TlD zWTTY (51) 



With these facts, expression (50) can be written as 

e#w-i e ftw-a ...gftv-j. (52) 



Z ( N 3) _ J m _i_R m _ 2 . . . M N _ 



Z'(N) 

In the thermodynamic limit iV->oowe can write e^^ -1 = . . . = e? N -i ~ e^, and therefore 

N 

y^j e -j/3(e k -it) 

(n k ) = ^ • (53) 

e-jffOfc-M) 

j=0 

can be readily evaluated to give the average number of particles as a function of temperature 
T and the chemical potential \i. The wanted relation reads 

^ = e ^-,)-r (54) 

which is the well-known Bose-Einstein (BE) distribution for an ideal gas of integer spin 
particles. In particular, the case considered here correspond to zero-spin particles. Note that 
taking the limit iV — > oo is a crucial step to obtain properly the Bose-Einstein distribution, 
since our starting point is the canonical partition function of exactly N particles. 



It is worth to notice that the chemical potential /x in this context arises as a consequence 
of a physical restriction: the constancy of the total number of particles. From that con- 
sideration it is not obvious or straightforward to see whether the chemical potential is a 
positive or negative quantity. In order (n^) be a non-negative quantity it is required that 
e i — A 4 — f° r a U 3- This implies that /j < €q, where eo is the single-particle ground-state 
energy. Since in general, eo — > in the thermodynamic limit, fi < for all temperatures. 
Note that the restriction imposed by Eq (45), gives an implicit definition of \l in terms of 
the particle density n = N/V and temperature T. This is, 

oo oo 1 

k=0 k=0 
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where the sum is strictly over an infinite number of energy levels since any particle can be 
occupy in principle any energy level accessible to the system. 



To obtain the equation of state /i = fi(n, T) implicitly defined in Eq (55) it is con- 
venient to transform the sum into a integral using the density of states p(e). Given the 
relationship between the wave vector A; of a free particle contained inside a box of volume 
L 3 (under periodic boundary conditions) given by quantum mechanics k x = 2irn x /L with 
n x = 0, ±1, ±2, . . ., the sum ^ fc over the wave vectors can be written in terms of an integral 
over the numbers n x ,n y and n z . A change of variable enable us to write that 

L V f „ V 



?^ujy*=^y* « 56) 

which can be represented in terms of the density of states p(e) as 

V [ p(e)de (57) 
k J 

where p(e) is a function that depends on both the system itself and its dimensionality 18 . In 
three dimensions, 

m 2/3 

with m the mass of the particles. 



For free bosons in three dimensions, Equation (55) can thus be written as an integral 
over the energy levels e as 

where h = h/2n and V the volume of the system. As usual (3 = (ksT)" 1 and e is the energy 
of the system. The integral in Eq. (59) can be expressed in terms of more familiar functions 
using a change of variable. For the sake of clarity we present this calculation in some detail. 
A new variable x can be defined clS X = /3e. Then, 

e^de 1 r x^dx 



l e^~^ - 1 f3 3 / 2 J (e x - 1 ^ 60) 
where ( is defined as ( = e _/3M . The last term in Eq (60) can readily identified with a 
particular type of special function. This is the Poly- logarithm function Li s (Q defined as 

r- 1 f°° tS ~ 1(lt 

^KHrwi C^T- (61) 
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For s = 3/2 and ( — > ( 1 we have 



, i\ i r°° t" 2 dt 

L ' 3 "{c) = nml —i (62) 



which enable us to rewrite the last integral in Eq (60) as 



00 x 



With this result the chemical potential /x can be written as an implicit function of temper- 
ature T and particle density n = N/V as 

m 3 /V/ 2 (A; B r) 3/2 



n 



-^3/2 (e^) . (64) 



2V2tt 2 ^ 3 

Equation (64) corresponds to the equation of state \i = /i(n, T) for an ideal Bose gas. 
This is completely equivalent to the standard equation of state for density n in terms of 
volume V and pressure p as has been shown in standard thermodynamics textbooks. Both 
/j = /i(n, T) and n = n(T,p) contain the same information and thus can be used indistinctly 
to obtain thermodynamic information of the system. In Fig. 3, the monotonic dependence 
on temperature of the chemical potential is shown for a fixed value of the density. 

As discussed in many textbooks 6,8 a phase transitions occurs at a critical temperature T c 
when fi = 0. From Eq. (64) that temperature is given by 

Tc = JlL-J—n 2 / 3 (65) 
C(3/2) k B m 1 ; 

marked with a dot in Fig. 3, its value in units of Tp is given by Tc/Tp = 
[4/(3C(3/2)2V2)] 2 / 3 ~ 0.436. In expression (65) the quantity ({3/ 2) = Li3/2(1) is the 
function zeta of Riemann. In Fig. 4, isothermal curves (light-color) of fj,(n,T) are shown. 
The critical density n c at which BEC occurs is determined by /i(n c ,T) = 0, and is given by 
n c = {mk B T/27ih 2 ) 2 ^C{3/2). 



The peculiarity that the single-particle ground state vanishes in the thermodynamic limit 
can be used to discuss a thermodynamic similarity between the photon gas and the uniform 
ideal Bose gas. Indeed, for temperatures smaller or equal to T c the particles in the ground 
state iVo (the condensate) do not contribute to the thermodynamics of the gas. Therefore 
we can disregard the condensate even when the total iV = N + N e is fixed, where N e 
denotes the number of particles occupying the single-particle excited states. Of course, the 
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FIG. 3: Chemical potential in units of Ef (the Fermi energy of a spinless ideal Fermi gas) as func- 
tion of temperature in units of Tp = Ef /ks, ks being the Boltzmann's constant, for: i) the ideal 
Fermi gas (continuous-blue line), ii) the ideal classical gas (red-dashed line), and Hi) the ideal Bose 
gas (magenta-dash-dotted line). The BEC critical temperature Tc/T F = [4/(3C(3/2)2\/2)] 2 / 3 ~ 
0.436 is marked with a dot while the Fermi Energy with a square. The inset shows how the chemical 
potential of the ideal quantum gas approaches the classical one at large temperatures. 

situation described only makes sense for the ideal case and is presented here just for academic 
purposes. Thus, the thermodynamic properties of the gas are dictated by the behavior of 
N e (T), which grows with temperature as occurs with the photon gas in a thermal cavity. 
The condensate plays the role of a particle source just as the walls of the cavity emits and 
absorbs photons from the cavity. The quantitative difference between both systems are 
the result form their different dispersion relations, = 4^- for the uniform Bose gas and 
ek = hck for the photons. To exemplify this, consider the temperature dependence of the 
specific heat at constant volume Cy, which for T < T c grows monotonically as T 3//2 for the 
uniform Bose gas and as T 3 for photons. This analogy would lead one to conclude that a 
photon gas behaves as a Bose gas with a infinite critical temperature. It is possible however 
to reduce the critical temperature of the photon gas to finite temperatures. One just have 
to manage to make the number of photons to be conserved at some critical temperature 
and this seems to have been recently realized experimentally by Klaers et al. by using an 
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ingeniously experimental setup 30 . 
B. Fermi-Dirac Distribution 

Regarding fermionic systems, the electron gas has been a paramount system in solid state 
physics since the crucial observation of the fermionic character of the electron. Indeed, the 
electron gas has played a fundamental role in the first stages of the theory of metals 32 and 
on the understanding of the stability of matter 33 . After some years of the experimental real- 
ization of condensation in a degenerate Bose gas 20 ' 21 , researchers started to turn their sight 
to the Fermi gas. The first experimental realization of a degenerate Fermi gas was carried 
out by de Marco and Jin 34 exhibiting the consequences of Pauli's exclusion principle. This 
work has triggered a renewed interest on Fermi systems, not only to the understanding of 
phenomena that emerge from strongly interacting fermion system in condensed matter, such 
as superconductivity, fermionic superfluidity etc., but also, to test and probe the theoretical 
predictions of quantum mechanics. In particular, the trapped ideal Fermi gas has been a 
system under intense theoretical research in the latest years 35-42 . Many other experiments 
have been developed to unveil the fermionic properties of matter 43-46 . 

For particles obeying Pauli's exclusion principle, the possible values of rij are restricted to 
and 1. In addition, a system of fermions, like the electron gas, must satisfy that total num- 
ber N of particles must be a constant, i.e., condition (45) must be satisfied. As mentioned 
before, restriction on the total number of particles implies that if a particular energy level 
is occupied by one particle, then the remaining N — 1 particles should distribute themselves 
into different energy levels. With this as the key idea we proceed to the calculation of the 
average number (n&) in the case of a Fermi-Dirac. We tackle this calculation in a different 
way from what we did in the Bose-Einstein case. For the FD we shall follow closer the 
procedure suggested in 7 . 

In a similar way as we did in the BE case let us define the sum 



n q ^n k 

as the sum for N particles carried out over all energy levels different from e&. The average 
number (n^) can thus be splitted into the /^-contribution and the remaining terms different 




(66) 
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from k. This is, 

Er,, p -0n k e h V p -^(nie 1 +n 2 e 2 + ...) 

\ k > V e -/3n fe € fc Y- e -/3(n iei +n 2 e 2 +...) " ^° ') 

/L^in k Z^n q ^n k 

For a system of fermions, can only have two values, zero or one. In addition since the 
number N of particles is fixed once a particular energy is populated with one particle, the 
remaining ones must be occupied by N — 1 particles, Eq(66) together with Eq (67) can be 
written for a fermion system as 

e fckZ'(N - 1) 

{Hk) = Z'(N) + eP*Z'(N-iy m 



In order to relate Z'(N) with Z'(N — 1) it is useful to consider the Taylor expansion of 
the quantity log Z'(N - AN) . For AN < N, 



If we define as 



log Z\N - AN) ~ log Z\N) AJV. (69) 



- dl °Z Z ' err* 
QN = ^N~' (70) 



we can write Eq (69) as 

log Z\N - AN) ~ log Z'(iV) - a N AN, (71) 

which yields to 

Z\N - AN) ^ 



Z'(N) 



e - aNAN . (72) 



Let us remember that Z'(N) is a sum defined over all states excepting the k one. One may 
expect then that for N ^> 1 variations in the logarithm may be some kind of insensitive to 
which particular state s has been omitted. Then, it may be valid that does not actually 
depends on the state k chosen and thus we can simply write = a 7 . Inserting this in 
Eq (72) and performing the sum in Eq (68) accordingly, we obtain for a fermi system the 
well-known Fermi-Dirac Distribution, 

< n *> = eftJa + V (73) 
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where a s is given formally by Eq (70). A direct interpretation for a can be given in terms 
of the chemical potential \i by recalling that 

Since F = —ksTlogZ, then a = —^jhsT- The Fermi-Dirac Distribution can be written 
then in a more usual form as 

< n *> = + V (75) 

where \x is up to this point an undetermined quantity that can be obtained by imposing the 
following condition 

s k 

Since both (rife) and N should be positive quantities the chemical potential \i must adjust 
its value in agreement with the value of the energy levels in such a way that N > and 
(nk) ^ be fulfilled in any physical situation. 

In the same spirit as we did for our calculation in the BE case, it is possible to go a step 
forward to calculate explicitly the equation of state \x = fi(n, T) for the ideal Fermi gas. 
By using expression (58) for the density of states, Eq. (76) can be written as the following 
integral 

N=^—(^) 3/2 r — (77) 
r(3/2) J e ^-M) + l' 

and in terms of the particle density n and the polylogarithm function Li s (z) we have 



n 



Li m {-e^), (78) 
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which gives, implicitly, the equation of state fi(n,T). In the zero temperature limit, the FD 
distribution (75) has a step-like shape 9(e — fi), where 6(x) is the Heaviside step function 
that takes the value 1 if x > and otherwise, thus, the chemical potential /x(n, T = 0) 
coincides with the so called Fermi energy Ep = ksTp = h 2 k F /2m whose dependence on n is 

E F = ^-j6n'nf S . (79) 

Due to the exclusion principle only one fermion can be allocated in a single-particle energy 
state (with no degeneration). Thus given iV particles, the system's ground state is obtained 
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by filling the first N single-particle energy states. The Fermi energy corresponds exactly 
to the last occupied state. For finite temperatures, but still much smaller than the Fermi 
temperature, the Fermi-Dirac distribution is modified from its zero temperature step-shape 
only around fi F ~ E F and the chemical potential can be computed by the use of the 
Sommerfeld approximation (see Ref. [47] for details) giving the well known result 

'2 



fi F 



E f 



l-'"-{T/T F ) 2 + ... 



(80) 



The temperature T* that separates the fi > region from the fi < one, can be computed 
exactly and is given by T* = [T(5/2)C(3/2)(1 - y/2/2)]~ 2 ^ T F ~ 0.989 T F , where T F denotes 
the Fermi temperature. 




(dark color) fermions 
(light color) bosons 
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FIG. 4: Chemical potential fi in units of ksT c for the ideal Bose and Fermi gases as function 
of the particle density n for various values of temperature. T c corresponds to the Bose-Einstein 
condensation critical temperature (65) of a boson gas with the arbitrary density no- Ao corresponds 
to the thermal wavelength evaluated at T = T C . Note that both cases converge to the same values 
of the chemical potential for small enough density, this corresponds to the classical limit. 



Observe that ft can be a positive quantity, even in the thermodynamic limit, in contrapo- 
sition to the Bose and classical gas where it is always a negative quantity. This behavior is a 
direct consequence of the quantum effects at low temperatures, in this case arises from the 
exclusion principle. At zero temperature, we can straightforwardly use the discrete version 
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of Eq. (1) to compute \i. This is so since AS = when adding exactly one particle to the 
system and therefore \i = AF = AU = Ep > 0. In Fig. 3, fi exhibits a monotonic decreasing 
dependence on temperature (blue-continuous line). Note that the transition to the classical 
behavior can occurs at very high temperatures, as high as the Fermi temperature which for 
a typical metal is of the order of 10 4 K. In Fig. 4 the dependence of /i on the particle density 
is shown for various isotherms, for this we have chosen the scaling quantities /j,q, hq and To 
of a reference system consisting of Nq particles in the volume V and Fermi energy given by 
(79). At low densities and finite temperatures [i is negative exhibiting the classical behavior. 

IV. QUANTUM STATISTICS BEYOND BOSE-EINSTEIN AND FERMI-DIRAC 

As reviewed in detail in previous sections, Bose-Einstein and Fermi-Dirac statistics de- 
scribe quantum systems of particles with complete different macroscopic thermodynamic 
effects. The essential difference between BE and FD systems is the Pauli exclusion principle 
which hinders the occupancy of a particular energy level to the values and 1. In view 
of this, we address the question: Is there any intermediate case between the BE and FD 
statistics? Recall that both can be viewed as extreme opposites of occupancy. Whereas 
BE enables any number of particles from zero to oo, FD blocks out any possibility beyond 
single-occupancy. In this Section we explore the possibility of Intermediate Quantum Statis- 
tics (IQS), i.e, statistics where any single-particle energy level can be occupied by at most j 
particles, with j an integer number between zero and oo. This is the most general case with 
BE and FD particular cases corresponding to j = oo and j = 1, respectively. 

Let us denote with IQS.,, the IQS of order j, of a non- interacting quantum system of 
particles where any single-particle energy level can be occupied, at most, by j particles. The 
calculation of the average number of particles {n^)j at the energy level k corresponding to 
the statistics IQS., can be done in a straightforward manner by generalizing the procedure 
used here to calculate the BE and FD ideal statistics. As an illustrative case let us consider 
the calculation of {nk}2 which is associated to a quantum system where the single-particle 
energy level can be occupied by zero, one or two particles. In general, as we have reviewed 
previously, the average number (rife) can be written as 




(81) 
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where the notation is exactly the same as before. The key issue to proceed with the 
calculation is to realize that the two sums on the numerator and denominator in Eq (81) 
are interrelated due to the restriction N = X]fc( n fc)- ^ or example, = 1 the sum 
Y2 n ^n k e _/3 * ni<El+n2€2+ ' ") must be performed over the N — 1 remaining particles since N is a 
fixed quantity. 



For a system obeying the IQS 2 statistics each energy level may be occupied by zero, 
one or two particles. Then, taking into account this (rik)2 can be written explicitly as 

_ e-^Z'jN - 1) + 2e-^Z { '\N - 2) 
^ k ' 2 ~ Z'{N) + e-faZ'{N - 1) + e~ 2 ^Z'{N - 2) ^ ' 

where 

Z\N) = e- p{ni€1+n2t2+ - ] (83) 

n q ^n k 

is a sum performed over iV particles leaving apart the energy level e&. Accordingly, Z'(N-l) 
represents the same sum performed over N — 1 particles, Z'(N — 2) a sum performed over 
N — 2 particles and so on. In general, Z'(N) and Z'(N — AN) are related at first order by 

Z'(N - AN) = Z\N)e~ aAN . (84) 

where a is the fugacity and is related to the chemical potential \i by a = — Then, (n k ) 2 
can be written as 

e -/3(efc-M) _|_ 2e~ 2l3 ( ek ~^ 

^ 2 = i + e -P{**-n) + e -Wte-n) (85) 

which represents the average occupancy for a quantum system with IQS2 statistics. It is 

worth to notice that the case IQSi, which represents the well-known FD statistics, is included 

in this expression. In such case the last terms in both the numerator and denominator are 

dropped out obtaining for IQSi 

e -P(e k -n) J 

fa** 1 = 1 + e -K*u-ti = e K* k -») + I ( 86 ) 
which is the Fermi-Dirac Statistics. 



The procedure outlined above can be readily generalized to calculate the average occu- 
pancy for a system with a IQSj statistics. In such case, the single-particle energy levels can 
be occupied by zero, one, two up to j particles simultaneously. This is the most general case 
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of a intermediate statistics between the Fermi-Dirac and Bose-Einstein cases. Please note 
that whereas the FD statistics corresponds to IQSi, the BE statistics lies on the opposite 
extreme where the occupancy j tends to infinity. The general expression for the average 
occupancy (n k )j in the IQSj case reads as 

(n ) = ^™-*toZ\N-r) 

£ J r=0 e-^Z'(iV-r) 1 ^ 

which can be calculated explicitly as 

e (l+j)(a+Pe k ) _|_ j _ M _|_ j\ e (a+/3e fc ) 
^ j = [ e (<*+Pe k ) _ 1] [ e (l+j)(a+0 £fc ) _ ( 88 ) 

where j can run from zero to infinity. In order to check out that this expression is correct 
let us calculate some particular cases. For j = we obtain the trivial limit case with no 
statistics at all. If no particles are allowed to occupy any energy level then there is no 
average occupancy. For j = 1 we recover the FD statistics since the expression for IQSi 
obtained directly from the substitution of j = 1 in Eq (88) 

e 2(a+/3e fe ) + -|_ _ 2 e ( a +P*k) 



is completely equivalent to Eq (86). The BE statistics can be also reproduced from Eq 
if we consider, as we did previously in the standard derivation of the Bose-Einstein statistics 
that the occupancy can run from zero to infinity. Then, the sums in Eq (87) must be carried 
out from zero to infinity. When that consideration is taken properly, (rik)j in the limit when 
j — > oo reproduces the BE case since 

where a as usual is the fugacity. 

The procedure outlined and described here to calculate the average occupancy (n^j in 
a quantum system obeying a IQS of order j is based on a procedure suggested by Reif 7 
for the calculation of BE, FD and Planck distributions exclusively. The generalization 
provided here shows that the chemical potential /i and its associated quantity the fugacity 
a are physical quantities related not only to the BE and FD statistics but to all types of 
statistics that preserve the total number of particles. The methodology proposed here can 
be straightforwardly explored with undergraduate and graduate students in order to clarify 
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how the concept of chemical potential arises and what is its role in the development of the 
standard FD and BE statistics. As an interesting issue to explore in this direction, it is 
worth to underline that once the restriction of the preservation of the number iV of particles 
is imposed this automatically restricts the summations implied in Z'(N), Z'(N — 1), . . .. All 
these sums are related and the connection factor is the fugacity a of the system. If these 
facts are not properly taken into account, all the sums defined by Z'(N), Z'(N — 1),... 
may be wrongly taken as the same. This misconception will bring the cancelation of the 
connection factor implied. 

To finalize this section we would like to make some comments on different approaches that 
have been proposed to deal with quasi-particles that are neither bosons or fermions. One of 
them is the concept of particles with fractional statistics also known as "anyons" introduced 
by Leinaas and Myrheim 48 and Wilczek 49 in two dimensional systems and that has found 
application in the theory of the fractional quantum Hall effect and anyon superconductivity. 
A completely new concept without reference to dimensionality was developed by Haldane 50 
based on the idea that the dimension D of the Hilbert space of single "particles" (in general 
quasi-particles that result from topological excitations in condensed matter) changes as 
particles are added to the system according to AD = —gAN. In other words, quantum 
correlations between "particles" are introduced by making the available states to depend on 
which states have been already occupied. The Bose statistics is recovered by setting g = 
and Fermi if g = 1 . 

V. WEAKLY INTERACTING QUANTUM GASES 

Let us finalize this brief review on the role of chemical potential in classical and quantum 
gases by briefly addressing the case of weakly interacting quantum gases. This case turns out 
to be of great relevance since it is the standard theoretical model to analyze Bose-Einstein 
Condensation in alkali atoms under magnetic and optical traps. The interacting Fermi gas , 
on the other hand, lies at the foundation of the superconductivity and fermionic-superfluidity 
theory when the effective interaction between fermions is attractive. 
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A. The Bose Gas 



In order to describe the dynamics of a weakly interacting Bose gas it is customary to 
start with the general Hamiltonian operator H given by 



H 



dr¥ ( - A- + V ext (r) | + 



J ' drd?ty\f^ty\P)U(r-P)ty(r)ifr(r') 



(91) 



where ty(f) and ^(r) are the field operators of annihilation and creation of particles at 
position r and U(f — r') is the interacting potential between two particles. In general, the 
experimental situations involve an external potential V ex t(r). In the case of bosons, the field 
operators satisfy a particular set of commutation rules given by 



S 3 (f{ - r 2 ) 



and 



0. 



(92) 



(93) 



The complete solution of Eq. (91) for any arbitrary potential U(f— r') is a formidable task 
beyond our current capabilities, however, for some particular situations it is possible to make 
a step further to approximate the potential U if — r 1 ) as a contact potential represented by 
a Dirac delta function 

U(f-r f ) = U 5 3 (r-r r ), (94) 

where Uq is the strength of the interaction given by Uq = 4ira s /m with a s the scattering 
length and m the mass of the particle. This has proved to be particularly accurate to 
describe interactions in Bose gases composed of alkali atoms like 23 Na, 87 Rb, 7 Li at very low 
densities and temperatures. In such systems, the interaction occurs via a s-wave quantum 
scattering process with a s the relevant parameter that characterizes the interaction between 
atoms. 

With these considerations it is possible to rewrite Eq (91) as 

U 



2V 



p>0,<?>0,r^0 

which is a second-quantization representation 52 in the momentum space q for the Hamilto- 
nian of the weakly interacting gas. The operators and a q are creation and annihilation 
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operator in the momentum space. The Hamiltonian in Eq. (95) can be split up into the zero 
momentum state q = and states with g ^ 0. Neglecting terms of the order N" 1 which 
vanish in the thermodynamic limit, Eq. (95) is written as, 

]\j2TT TT 

F = ^ + E ( e S + 2n ° u °) h % + y E a&afrr. (96) 

q^O p,q 

where no = Nq^q 2 is the density of the condensate and \l/o is the corresponding wave func- 
tion. This Hamiltonian can be expanded around an equilibrium occupation distribution f q 
which for bosons is the Bose-Einstein distribution function 5 . To first order, the Hamiltonian 
is 

9^0 p,q 

where n = n + n\ is the total particle density of the system and n x = J2 q >o ^q^q * s the 
density of uncondensed particles which is a sum over all the non zero momentum states. The 
Hamiltonian in Eq (97) is known as the Hartree-Fock (HF) approximation for the weakly 
interacting Bose gas. The second term of Eq (97) shows the intrinsic nature of the Hartree- 
Fock approximation as a mean field theory. The energy e q + 2nUo to add or remove a particle 
to a state with non zero momentum is an average over all the pairwise interactions between 
particles. 

The equations for the wave functions of the condensate ^o(r) and the uncondensed phase 
^^(r) can be obtained from the Heisenberg Equation —i/h[H,^f] = d t ^ with H given by 
Hartree-Fock approximation in Eq (97) as 

ti 2 

-— V 2 * + 2m Wo + n U ^ + V ext ^ = eo^o 
2m h i (98) 



-— V 2 * q + 2U [m + n ] + V ext V q = e q ty 



2m 



q 



with the last equation valid for q ^ 0. Since N and Nj. for k > are assumed to obey 
a Bose-Einstein statistics, Eqs. (98) enable us to obtain the thermodynamic framework of 
the interacting Bose gas in the HF approximation. These equations determine the chemical 
potential \x in terms of the total particle density n and the temperature T as 

n\ 3 T = g 3/2 [P(v-2nU )] (99) 

for T > T c and 

n = n + i# 3 /2 (-pnoUo) (100) 
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with n = C/ (2n — no) for T < T c , with 



9v{a) 



T(v) 



- 1 



-dx 



(101) 



is the Bose integral, and At is the thermal de Broglie wavelength (31). For a system of units 
where a s = m = h = kg = 1, Eqs. (99) can be written in dimensionless form as, 

3/2 



n 



T 
2^ 



#3/2 



T 



for T > T c and 



rp \ 3/2 

" = "» ( ^ ) ^/2 



47m 



(102) 



(103) 



with /i = 47r (2n — no) for T < T c . Here the strength of the interaction Uq has been replaced 
by its dimensionless An value. 



Eqs. (102) and (103) contain all the relevant thermodynamic information for the weakly 
interacting Bose gas in the HF approximation and it is possible to solve them for \x in terms 
of the total particle density n and temperature T in order to obtain the isotherms of the 
equation of state fi = /i(n, T) for a gas confined in a box of volume V. Recently one of us has 
address this issue 53 for different values of the gas parameter 7 = a s 3 n obtaining isotherms 
for the weakly interacting Bose gas in the HF approximation (Fig. 5). The results show 
that the HF approximation while a valid theory of the interacting gas near zero temperature 
fails to predict and adequate physical behavior near the transition. Indeed, in the vicinity 
of the critical density n c , the HF formalism predicts a non single-valued profile for w(n) a 
feature forbidden by fundamental thermodynamic principles. 



B. Interacting Fermi gas 

The corresponding Hamiltonian for fermions with two spin states a =7, J,, as naturally 
occurs in several condensed matter systems, may be written as 

H = J2[ rfr ^W (-^ + V ext (a,r)-A * CT (r) 

+ y*drdi^*{(r)*J(r / )?7(r-r')*4.(r)* t (r / ) (104) 



33 




FIG. 5: Chemical potential \i as function of total density n in the Hartree-Fock approximation 
for two fixed temperatures ksT = 0.01 and k^T = 0.1 in units where m = H = a s = 1. For 
temperatures T above the critical temperature T c , i.e, densities n below the critical density n c 
the HF approximation gives correctly the behavior of the chemical potential for which [i —> — oo. 
However, for T < T c or equivalently for n > n c , the HF approximation yields to a chemical 
potential with a non single-value behavior. This is unacceptable based on fundamental principles 
of thermodynamics. 

where the field operators obey the fermionic anticommutation relations {^(r), ^,(r')} = 
^o.(r)^^,(r') + ^^. / (r / )^' -(r) = 5 a y5(r — r'). The density of fermions in each spin state 
n-fg, fixes the energy scale Ep of the noninteracting fermion given by expression (79) with 

n = i n t = \ n i- 

Two different aspects of the interacting Fermi gas are obtained depending on whether 
the interaction between fermions is attractive or repulsive. At zero temperature cor- 
rections to the ideal case value of the chemical potential can be calculated in terms 
of the scattering length a s which measures the interaction strength. In the dilute limit 
we have /ci?a s 1, where the Fermi wavevector /cp depends on the particle density as in (79). 
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In the repulsive case there is no change in the intuition we have developed from the 
interacting classical gas. So one would expect the chemical potential to rise above the value 
of the noninteracting case. This is clear from the Landau theory of the Fermi liquid 54 . A 
calculation due to Galitskii (see Ref. [52] pp. 147) gives that for kpa sufficiently small, the 
chemical potential can be written as 



exhibiting that repulsive interaction leads to an increase of the chemical potential as it occurs 
in the classical case. There is a particular interest in the case when attractive interactions 
between fermions of different spin-polarization are considered. This is, for example, the 
case in superconductors or in ultracold Fermi gases manipulated through magnetic fields, 
where the possibility of Cooper pairing is present. The formation of Cooper pairs, even for a 
extremely weak attraction, makes the Fermi sea unstable 58 leading to a ground state different 
from the Fermi liquid one called the Bardeen-Cooper-Schrieffer (BCS) ground-state. This 
microscopic mechanism developed further by Bardeen, Cooper and Schrieffer 59 serves as the 
basis for the explanation of conventional superconductivity. 

In the weak coupling limit, the chemical potential at zero temperature does not deviate 
significantly from the the noninteracting case value Ep. This picture changes if the strength 
of the attractive potential is increased and \i is computed in a self-consistent theory called 
the BCS-BEC crossover. Such a theory, introduced by Eagles 60 and Leggett 61 , extends the 
BCS one, where now the chemical potential changes due to the formation of tightly bound 
fermion pairs. The BCS-BEC crossover has been developed and widely applied, first in 
the context of high-Tc superconductivity and most recently in the formation of fermionic 
molecules in ultracold Fermi gases. We restrict our discussion to the case of zero temperature 
giving references for the finite temperature case. 

The Hamiltonian (104) can be written in momentum space as 



where Vfc,fc' is the two-body attractive interaction. In a self-consistent mean-field theory, 




(105) 



H 




(106) 




called the "gap" obeys the well known gap 



equation 




(107) 
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where 




[(e k -/,) 2 + A^ 



Ok - ^) 



(108) 



gives the quasi-particle density with wavevector k. The simultaneous solution of these 
equations give fi and A as function of the strength of the interfermionic interaction. 

The BEC-BCS crossover has implications on the behavior of fi as a function of the 
interaction strength 56 . In the weak coupling regime, /i = Ep and ordinary BCS theory 
applies. At sufficiently strong coupling, \x starts a monotonous decreasing behavior with 
increasing the attraction strength, eventually it crosses the zero value and then becomes 
negative in the Bose regime. With the appearance of tightly-bound pairs but // still positive, 
the system has a remnant of the Fermi surface, and we say that the many-body system 
preserves a fermionic character. For negative /i, however, no trace of a Fermi surface is left 
and the system is considered rather bosonic. 

VI. CONCLUSIONS AND FINAL REMARKS 

The concept of chemical potential in the context of classical and quantum gases has 
been revisited. For the classical gas, we started on giving a physical argument on why one 
should expect fi < and then we considered the effects of the inter-particle interactions for 
the particular case of the van der Waals gas. Based on a equivalence due to Widom that 
relates the work necessary to bring an extra particle from infinity to a given position in the 
system, with the difference in the chemical potential respect to the perfect gas, we give a 
physical argument on what to expect when interactions are considered. Thus, in the case 
of a pure repulsive interaction of the hard-sphere type of radius d, the chemical potential is 
shifted above the ideal gas value by the amount ^7i(d/l) 3 kBT. This result is valid only in 
the low density limit, i.e., d <^ I. For this case it is intuitively clear, that in the situation 
just described, it must be spent more energy in bringing an extra particle to the system 
than when no interactions are present at all and that this amount of energy increases as 
the density does. Thus an infinite amount of energy will be required to add a particle 
to a high dense classical gas. As the reader can now expect, attractive interactions shifts 
the chemical potential below the value of the non-interacting case. This is the case when 
considering the attractive tail in the model potential (23). 
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For the ideal quantum gases, we present a pedagogical way to obtain the Bose-Einstein 
and Fermi-Dirac distributions starting from a canonical-ensemble calculation. In the case 
of bosons, we discussed the implications on the chemical potential when considering a 
system, both, with conserving and non- conserving number of particles. In the later case, 
no reference to \x is needed, however, /i appears naturally once the number of bosons is 
required to be conserved. In this case n(T) decreases monotonically with temperature lying 
below, but asymptotically approaching, to the classical curve (see Fig. 3). In contrast, 
for fermions it was shown that /x acquire positive values due to the statistical correlations 
induced by Pauli's exclusion principle. In addition, fi(T) decreases monotonically from the 
Fermi energy lying above, but asymptotically approaching, to the classical curve (see Fig. 
3). On the light of these observations, we can use the ideas exposed for the interacting 
classical gas. Indeed, if we consider the quantum gas as classical, with quantum correlations 
given by a statistical interparticle potential v s t a t,ij 6 ' 62 , then, due to the attractive/repulsive 
nature of v s t a t,ij for the Bose/Fermi gas, their respective chemical potentials vary with 
temperature below/above the classical one. 

We have also briefly discussed the consequences of considering an extension of the ex- 
clusion principle when a single-particle energy level can be occupied at most for j particles. 
Finally, we have presented a discussion on the behavior of n for the case of weakly inter- 
acting quantum gases. For the Bose gas, fj,(T,n) gives important information on the nature 
of the BEC phase transition. In the case of the attractively interacting Fermi gas, fi gives 
important information on the nature of the system as the interaction strength is varied, 
going from loosely bound pairs (Cooper pairs) in the weak coupling to bosonic thightly 
bound-pairs in the strong interaction limit. 

VII. APPENDIX 

The basic idea to evaluate the configurational integral Qn given by Eq. (17), is to com- 
pute the statistical average of e~P Vl > 2 e~P Vl > 3 ■■■e~" VN - 1 < N over all possible configurations of 
the particles positions denoted with Qn/V n = e~P Vl > 2 e~P Vl > 3 ■ ■ ■ e~P v * f - l > N . van Kampen's ap- 
proach is based on a factorization of D into terms that takes into account the correlations 
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of k > 2 particles, i.e., 



N 



Qn/v n = JJ(4)(*) , (109) 

fc=2 

where (^) gives the number of combinations of k particles taken from the total N, and 

(110) 



-Pvi,2p—Pvi,3 . 



D 



with D is the immediate lower approximation for the same numerator. 

For k = 2, g?2 = e~P Vl ' 2 since D = 1 in this case. Thus the first factor in Eq. (109) is 
given by 

r dr 2 e~^ 2 . (Ill) 



N{N-l)/2 



V- 1 j d^V- 1 

In order to take the thermodynamic limit N, V — > oo with N/V constant, consider the 
following identity 

(N-l)/2 



V 



-/3ui,2 



V 



\ 1 N I ■ < 

1+ Nvl dT (e 



•j8«(r) 



1) 



(JV-l)/2 



;ii2) 



thus giving as result d\ 2 = exp < f^-Bi f with B\ = J dr (e -|8 "( r ) — l) . For dilute enough 
systems where only correlations of two particles are important this approximation should 
work fine. 

The calculation of the general factor (dk)^ k > is more involved and we present only a sketch 
of it. By writing = 1 + f id Eq. (110) can be rewritten as 



dk 



1 + /l,2 + ■ ■ ■ + f 1,2 f i,3 ■ ■ ■ fk-l,k 
D 



(113) 



van Kampen argues that the class of terms in the numerator of (114) that involve less than 
k particles and those that involve k particles but are reducible, are also present in D, such 
that the numerator can be written as (1 + X]{fc} A, 2/1,3 ' — l~ 0(V~ k ))D, where the summa- 
tion extends over all irreducible terms that involve 2, . . . , k particles. A term of the form 
I ' ' ' I Yli<j 9i,j d r i ' ' ' dr^, with g^j an arbitrary function of |r$ — rj | , is said to be irreducible 
if cannot be factorized into products of integrals of g^j involving less than k particle-position 
coordinates. For instance, it is straightforward to check that JJJ drxdr^dr^ gi t 29i,392,3 is ir- 
reducible while the integral JJJ dridr 2 dr 3 gi,29\,3 — V \J <irg(r)] 2 is not, here we have used 
the identity J dv — V. 
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Thus, we have 

d k = l+ { -^-^B k _ 1 + 0(V- k ), (114) 

where B k is given by (19) and we have recognized /i,2/i,3 • • • with the usual irreducible 
cluster integral 6-8 yF^ -Bfc-i- In the thermodynamic limit the factor (c4)( fe ) can then be 
written as 

and by combining this result with the result for k = 2 we finally get the desired result given 
by expression (18). 
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